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ABSTRACT 

Context. Studies of high-mass protostellar objects reveal important information regarding the formation process of massive stars. 
Aims. We study the physical conditions in the dense core and molecular outflow associated with the high-mass protostellar candidate 
IRAS 18566+0408 at high angular resolution. 

Methods. We performed interferometric observations in the NH 3 (J,K)=(1,1), (2,2) and (3,3) inversion transitions, the SiO J=2-l and 
HCN J= 1-0 lines, and the 43 and 87 GHz continuum emission using the VLA and OVRO. 

Results. The 87 GHz continuum emission reveals two continuum peaks MM-1 and MM-2 along a molecular ridge. The dominant 
peak MM-1 coincides with a compact emission feature at 43 GHz, and arises mostly from the dust emission. For dust emissivity 
index /? of 1.3, the masses in the dust peaks amount to 70 M Q for MM-1, and 27 M Q for MM-2. Assuming internal heating, the central 
luminosities of MM-1 and MM-2 are 6 x 10 4 and 8 x 10 3 L Q , respectively. 

The SiO emission reveals a well collimated outflow emanating from MM-1. The jet-like outflow is also detected in NH 3 at velocities 
similar to the SiO emission. The outflow, with a mass of 27 M Q , causes significant heating in the gas to temperatures of 70 K, much 
higher than the temperature of < 15 Kin the extended core. Compact (< 3") and narrow line (< 1.5 km s~') NH 3 (3,3) emission features 
are found associated with the outflow. They likely arise from weak population inversion in NH 3 similar to the maser emission. 
Toward MM-1, there is a compact NH 3 structure with a linewidth that increases from 5.5 km s -1 FHWM measured at 3" resolution 
to 8.7 km s -1 measured at 1" resolution. This linewidth is much larger than the FWHM of < 2 km s -1 in the entire core, and does not 
appear to originate from the outflow. This large linewidth may arise from rotation/infall, or relative motions of unresolved protostellar 
cores. 

Key words. ISM: kinematics and dynamics — ISM: H II regions — ISM: clouds — Masers — Outflows — ISM: individual 
(IRAS 18566+0408) — stars: formation 



•<-j . 1. Introduction IRAS 05358+3543, Beuther et al. 2002d; AFGL 5142, Zhang et 

' al. 2007). In the meantime, high resolution images in NH3 and 

»H ; Systematic surveys in the past decade identified hundreds other dense moleC ular gas tracers reveal interesting kinematics 

of high-mass protostellar candidates (Mohnan et al. 1996; close to massive protostars (Zhang et al. 1998, 2002). 
Sridharan 2002; Fontani et al. 2005). These objects, selected 

initially from the IRAS point source catalog, typically have In this P a P er < we P resent a hi 8 h resolution study toward the 

far infrared luminosities of >10 3 L , contain 10 2 - 10 4 M of high-mass protostellar candidate IRAS 18566+0408. At a kine- 

dense molecular gas (Molinari et al. 2002; Beuther et al. 2002a; matlc distance of 6.7 kpc (Sridharan et al. 2002), the source has a 

Williams, Fuller & Sridharan 2004; Beltran et al. 2006), and are far infrared luminosity of several 10 4 L . The object was initially 

associated with massive molecular outflows (Zhang et al. 2001, undetected at 2 and 6 cm at an rms of 0.16 mJy and 0.1 mJy, 

2005; Beuther et al. 2002b). Compared with ultra compact Hn respectively (Miralles, Rodriguez & Scahse 1994), but later de- 

(UCHn) regions, high-mass protostellar candidates have similar tected at 3 6 cm at a flux den sity of 0.7 mJy (Carral et al. 1999), 

amounts of dense molecular gas, but are less luminous and have and at 2cm at a flux of °- 7 mJy (Araya et al. 2005). 
much weaker emission at centimeter wavelengths. Therefore, This region is associated with H 2 maser emission at 22 

they are likely to be in an earlier evolutionary stage than the GHz, CH3OH maser emission at 6.7 GHz, and H 2 CO maser 

UCHn phase. emission at 8 GHz (Miralles, Rodriguez & Scalise 1994; Slysh 

These surveys were carried out mostly using single dish tele- et al. 1999; Beuther et al. 2002c; Araya et al. 2005). Dense 

scopes with angular resolutions of > 10". High angular res- gas traced by CS and CH3CN, as well as (sub)mm continuum 

olution imaging is required to probe dense cores and molec- emission is observed (Bronfman et al. 1996; Sridharan et al. 

ular outflows at spatial scales relevant to massive protostars. 2002; Beuther et al. 2002a; Williams et al. 2004). NH3 (1,1) 

In the past few years, images from (sub)mm interferometers and (2,2) emission was detected first by Miralles, Rodriguez & 

have often resolved poorly-collimated outflows identified by sin- Scalise (1994), and also by Molinari et al. (1996) (source 83) 

gle dish telescopes into multiple well-collimated outflows (e.g. and Sridharan et al. (2002) with single dish telescopes. 
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Beuther et al. (2002b) report a CO outflow in the northwest- 
southeast direction. The geometric center of the outflow, how- 
ever, is about 10" north of the 1.2 mm emission peak. Since the 
1 .2 mm continuum position is consistent with that of the submm 
emission (Williams et al. 2004), this offset is possibly caused 
by pointing problems in the CO observations with the IRAM 
30 telescope (Beuther, H., private communication). The 1.2 mm 
emission shows an extension of 15" toward the northwest of the 
peak emission. SiO J=2-l emission is detected in the region with 
a linewidth of 30 km s _1 at zero intensity (Beuther, H., private 
communication). 

The high resolution observations with the VLA and OVRO 
in this paper reveal a collimated outflow in SiO and kinematics 
in the dense core. In Section 2, we describe details of observa- 
tions. In Section 3, we present the main observational results. In 
Section 4, we discuss the different kinematic components in the 
region. A summary is given in Section 5. 

2. Observations 

2.1. VLA 

2.1.1. NH 3 observations 

The VLAQ observations of IRAS 18566+0408 were first con- 
ducted on 2001 July 23 in the NH 3 (J,K)=(1,1) and (2,2) 
lines in the C configuration. To improve the S/N in the data, 
follow-up observations were made from 2001 October to 2003 
January, in both CnB and DnC configurations in the NH3 (1,1), 
(2,2) and (3,3) inversion transitions. The integration time on 
source was typically less than 1 hour for each line. The point- 
ing center of the observations was RA (2000) = 18' ! 59'"09 5 .88 
and DEC(2000)=+04°12'13".6. We used 1849+005, 3C286 and 
3C273 as the gain, flux, and bandpass calibrators. The detailed 
parameters of the observations are summarized in Table 1 . 

The visibility data were calibrated using the NRAO 
Astronomical Image Processing System (AIPS). The uncertainty 
in the flux calibration is about 10%. The calibrated visibilities 
from different epochs were combined for the same line and im- 
aged in MIRIAD. The rms noise in the (1,1), (2,2) and the (3,3) 
lines is about 2 mJy in a 3" to 4" synthesized beam per 0.6 
km s wide channel. 

2.1 .2. 23 GHz and 43 GHz continuum 

The continuum observations at 23 GHz and 43 GHz were carried 
out with the VLA on 2002 September 26 and 2003 February 04, 
respectively. At 43 GHz, we used the fast switching calibration 
scheme that alternated between IRAS 18566+0408 and the gain 
calibrator 1849+005 in a cycle of 2mins. The total on-source 
time for IRAS 18566+0408 was about 2 hours at 43 GHz, and 
1 hour at 23 GHz. Calibration and imaging were performed in 
AIPS. The flux calibration was done by comparing to 3C286. 
The absolute flux scales are accurate to about 10%. The rms is 
0.1 mJy in the 43 GHz image, and 0.14 mJy in the 23 GHz im- 
age, respectively. 

2.2. OVRO 

The OVRO observations of IRAS 18566+0408 were carried out 
during 2002 November to December. The SiO J=2-l (v=0), 

1 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with the 
National Science Foundation. 



HCN J=l-0 and HCO + J=1-0 lines were observed simultane- 
ously in the lower sideband, along with 87 GHz continuum. The 
SiO line was observed with a total bandwidth of 3 1 MHz and a 
spectral resolution of 0.5 MHz (1.7 km s~'). The HCN line was 
observed with a bandwidth of 30 MHz at a resolution of 1 MHz 
(3.4 km s~'). The HCO + line was observed with a bandwidth 
of 7.5 MHz at a resolution of 0.25 MHz (0.8 km s 1 ). In ad- 
dition, the analog correlator provided continuum measurements 
of 4 GHz bandwidth. The pointing center of the OVRO observa- 
tions was the same as that of the VLA observations. The detailed 
parameters of the observations are summarized in Table 1 . 

The visibility data were calibrated in the OVRO MMA pack- 
age and exported to MIRIAD for imaging. The rms is 17 mJy 
per 1.7 km s _1 channel for the line images, and 0.5 mJy for the 
continuum. The HCO + emission is extended and suffers missing 
short spacing fluxes, thus, is not presented in this paper. 

3. Results 

3.1. Continuum Emission 

No emission is detected at 23 GHz at a lcr rms of 0.14 mJy. 
Figure 1 presents images of continuum emission at 43 GHz 
(or 7.0 mm) and 87 GHz (or 3.4 mm). The 43 GHz emission 
shows a compact feature with the peak position at RA (2000) 
= lS'^'W^, Dec (2000) = 04°12'15".7. The emission ap- 
pears to be slightly resolved with the 2".7 x 1".3 beam at a 
position angle of -7.6°. The peak and integrated flux densi- 
ties are 1 mJy/beam and 1.7 mJy, respectively, with a lcr er- 
ror of 0. 1 mJy. The emission has an extension in the northeast- 
southwest direction, which appears to be different from the po- 
sition angle of the beam. The deconvolved size of the emission 
is 2".0 x 1".2 with a position angle of 24°. This corresponds to 
a size along the major axis of 1 .3 x 10 4 AU. 

The 87 GHz emission is resolved with a 4". 9 x 4". 5 beam 
using natural weighting. The emission consists of a dominant 
peak, MM-1, coincident with the peak of the 43GHz emission 
to better than 0".l, and a secondary peak, MM-2, at RA (2000) 
= 18 /! 59'"09- , .21, Dec (2000) = 04°12'22".6. The peak flux den- 
sity of MM- 1 is 1 8 mJy/beam, with an integrated flux density of 
31 mJy. MM-2 is much weaker, with a peak flux density of 2.6 
mJy/beam. The uncertainty in these measurements is about 15%. 
There appears to an extented filament at a position angle of -58° 
connecting MM-1 and MM-2, which is better seen in the lower 
resolution (11") 1.2mm map in Beuther et al. (2002a). 

Figure 2 shows the spectral energy distribution of the con- 
tinuum peak MM-1. The 1.3 cm 3<x upper limit (lcr = 0.14 
mJy) is from this paper. The 3.6 cm continuum detection is from 
Carral et al. (1999). The 2 cm and 6 cm data are from Miralles 
et al. (1994) and Araya et al. (2005). The 1.2 mm, 850 fj.m and 
450 /im measurements are from the IRAM 30-m telescope and 
JCMT (Beuther et al. 2002c; Williams et al. 2004). The mid to 
far infrared data are from IRAS, MSX and Spitzer IRAC mea- 
surements. 

For emission at longer cm wavelengths, the contribution 
from dust is negligible. Toward MM-1, a faint continuum source 
was detected at 2cm (0.7 mJy) by Araya et al. (2005), and at 3.6 
cm (0.7 mJy) by Carral et al. (1999). However, we fail to detect 
the source at 1.3cm at an angular resolution of 1" and a lcr rms 
of 0. 14 mJy. The detections by Araya et al. (2005) and Carral et 
al. (1999) were made at resolutions of a few arcseconds, and can 
be reconciled with the non detection in this paper if the source 
is more extended than 1". However, an inconsistency remains at 
2cm at which Miralles et al. (1994) failed to detect the source 
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with a lcr rms of 0. 16 mJy at a resolution of 5". A possible rec- 
onciliation is that the flux varies with time. Despite the apparent 
differences, the faintness of the cm emission indicates that the 
massive star in this region is still extremely young in its evolu- 
tion, and has not produced significant free-free emission. 

The measurements from wavelengths shortward of 1.2 mm 
have poorer spatial resolution and sample a much larger area in 
the region. We fit a greybody model to the entire spectral en- 
ergy distribution from radio to infrared (IR) wavelengths. The 
far-IR measurements have a typical resolution of ~ V . The over- 
all spectral energy distribution can be fitted by three dust compo- 
nents, with temperatures of 210, 58 and 30K, respectively. The 
total luminosity of the region determined mainly by mid to far 
IR data at ~ 1' resolution is 8 x 10 4 L Q . The fluxes at the mm 
and submm wavelengths give a spectral index a of 3.9, defined 
asF v oc V ff , or/?= 1.9. 

For the compact continuum source MM-1, we use the high 
resolution 7 mm and 3 mm data to derive a power law index more 
appropriate for a mass estimate. To minimize the difference in 
beam size between the two frequencies, we image the 87 GHz 
data with a uniform weighting of the visibilities and obtain a 
peak flux density of 10 mJy/beam with a 3".3x2".0beam. These 
two values (10 and 1.0 mJy/beam) produce an upper limit to the 
spectral index a of 3.3, or an upper limit to the emissivity index 
P of 1.3. Assuming that the dust reaches an equilibium with the 
gas through collision at this high density environment (Burke & 
Hollenbach 1983), we approximate the dust temperature by the 
gas kinetic temperature of 80 K measured in NH3 (see Sections 
3.2, 3.3 and 4.1). For a dust opacity law k oc V 3 , and *-(250yum) = 
12 cm 2 g _1 from Hildebrand (1983), we obtain a mass within the 
4" .9 x 4" .5 beam (~ 30,000 AU) of 70 M G for (3 = 1.3. This 
mass is a small fraction of the mass (2 x 10 3 M ) estimated from 
the 1.2mm emission for the entire region (Beuther et al. 2002a). 

For the continuum peak MM-2, the non detection at 7mm 
gives a 3<x upper limit of 0.3 mJy. This value and the peak flux 
density of 2.6 mJy/beam at 87 GHz yield a spectral index of >3. 
Using assumptions similar to those for MM-1, we estimate the 
mass in the MM-2 core. With a temperature of 30K derived from 
the NH3 emission, we obtain a mass of 27 M for (3 = 1.3, and 
13M for/?= 1. 

3.2. Line Emission 

Figure 3 presents the integrated emission of the NH3 
(J,K)=(1,1), (2,2), (3,3) lines obtained from the VLA, and the 
SiO J=2-l and HCN J=l-0 transitions obtained from OVRO. In 
the NH3 (1,1) and (2,2) lines, there appears to be extended emis- 
sion in the northwest-southeast direction over a scale of 40". 
This component is relatively cold as the (2,2) emission is less 
extended than the (1,1) line. MM-2 is associated with NH3 gas 
and lies in a molecular ridge connecting MM-1 and MM-2 that 
is also seen in the dust emission. In addition to the extended 
gas component, a compact NH3 emission component associated 
with MM-1 is seen in all three NH3 lines, with a deconvolved 
size of 1".2. The peak intensity of the integrated emission in the 
(2,2) line is 0.22 Jy km s~ '/beam, similar to the value in the ( 1 , 1 ) 
line. Thus the gas in this component is relatively warm. In the 
NH3 (3,3) line, there is also extended emission in the northwest- 
southeast orientation. The extension, at a position angle of 135°, 
is similar to the emission in SiO and HCN. Unlike in the case of 
the NH 3 (1,1) and (2,2) emission, MM-2 is located toward the 
edge of the NH 3 (3,3), SiO and HCN emission. 

To show detailed kinematics in the region, Figure 4 presents 
the channel maps of the NH3 (J,K)=(1,1), (2,2) and (3,3) lines. 



The extended emission is present mostly at velocities from 84 to 
87 km s _1 , with a peak velocity of 85.2 km s _1 corresponding 
to the cloud systemic velocity (Bronfman et al. 1996). The typi- 
cal linewidth for the extended emission is about 1-2 km s _1 in 
FWHM and the typical temperature is < 15K (see the tempera- 
ture map in Figure 6 and discussions in Section 4). The relatively 
narrow linewidth and low temperature in the gas indicate that 
this extended component is from the quiescent gas in the core. 

At velocities of 86 to 87 km s _1 , there appears to be a molec- 
ular ridge (position angle of 148°) between MM-1 and MM-2 in 
the (1,1) and (2,2) emission. MM-2 coincides with a peak in the 
NH3 emission (see channel 86.4 km s~'). At velocities less than 
84 km s _1 and greater than 87 km s -1 , there appears to be com- 
pact NH3 emission toward the position of MM-1. This compact 
emission, with a FWHM of 5.5 km s -1 measured at 3" resolu- 
tion, is strong in the (2,2) emission relative to the (1,1) emission, 
indicating that the NH3 gas is rather warm. The ratio of the NH3 
(1,1) and (2,2) lines gives a rotational temperature of 45 K. 

In the NH3 (3,3) line, the extended component seen in the 
(1,1) and (2,2) is not as dominant. A compact emission com- 
ponent toward the position of MM-1 stands out prominently. 
Since the (3,3) transition has a higher upper energy level (124K) 
as compared to the (1,1) (23K) and (2,2) (65K) lines, the (3,3) 
emission confirms that the compact component is rather warm. 

In addition to the compact component toward MM-1, there 
appear to be four additional compact emission components in 
the NH3 (3,3) line, two toward the east of MM-1 in the velocity 
channels of 84.6 and 85.5 km s _1 , one to the west of MM-1 in 
the velocity channel of 85.2 km s _1 , and one to the northwest of 
MM-1 from velocities of 85.5 to 85.8 km s _1 . We refer to these 
features as 'A', 'B', 'C and 'D', respectively (The crosses in the 
channel maps in Figure 4). There appear to be no corresponding 
emission peaks in the NH3 (1,1) and (2,2) lines. Unlike the broad 
NH3 (3,3) line emission toward MM-1, these four components 
have rather narrow velocity width of about 1 .5 km s -1 in FWHM, 
but extended line wing emission. We will discuss these features 
further in Section 4. 

Besides the compact emission components in Figure 4c, the 
NH3 (3,3) emission also shows an extended structure in the 
southeast-northwest direction. This structure has rather broad 
line wings (15 km s _1 from the cloud velocity), and high tem- 
peratures of 70 K. MM-2 is not associated with any peaks of the 
(3,3) emission. 

Figure 5a presents the channel maps of the SiO J=2-l transi- 
tion. The SiO emission is elongated and lies mostly to the north- 
west of MM-1 at a position angle of 135°, similar to the extended 
emission in the NH3 (3,3) line. The SiO emission is present from 
velocities of 70 to 93 km s _1 . There appears to be higher velocity 
SiO emission toward MM-1, but none toward MM-2. Since SiO 
abundance is typically low in quiescent clouds (Ziurys, Griberg 
& Irvine 1989) and is enhanced by a few orders of magnitude in 
outflows (e.g. Zhang et al. 1995) due to shock processes (Pineau 
des Forets, Flower & Chieze 1997), the SiO emission here most 
likely traces a well collimated outflow originated from MM- 1 . 
The orientation of the SiO outflow is consistent with the bipolar 
CO outflow reported by Beuther et al. (2002). 

Figure 5b presents channel maps of the HCN 1-0 emission. 
The HCN 1-0 transition has three hyperfine components (F=l- 
1, 2-1 and 0-1) at relative frequencies corresponding to 4.8, 
and -7.1 km s _1 , respectively. We set the hyperfine component 
F=2-l at a Vlsr of 85.2 km s _1 , the cloud systemic velocity. It 
appears that the emission from the F=l-1 component is weak. 
The F=2-l and 0-1 components are detected around 80 and 73 
km s _1 , respectively (see also Figure 7), 5 km s _1 blue shifted 
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from the cloud velocity. The HCN emission arises mainly in two 
strong peaks: One is associated with the dust peak MM-1, the 
other is 4" offset from MM-2 and coincides with the SiO peak 
in the outflow. Little emission is detected toward MM-2. 

The measured flux ratios of the three HCN hyperfine com- 
ponents (F=0-1, 2-1 and 1-1) amount to 1:1:0.2. Under the LTE 
condition, line ratios vary from 1:5:3 for optically thin emis- 
sion to 1:1:1 for optically thick emission. The measured ratios 
are not consistent with partially optically thick gas under LTE. 
Anomalous hyperfine ratios of HCN have been found in both 
cold dark clouds and warm clouds around Hn regions (Walmsley 
etal. 1982;Cernicharoetal. 1984). Possible causes involveover- 
lapping hyperfine components of higher rotational transitions or 
core-envelope density structures/velocity gradients in the cloud 
(Gonzalez-Alfonso & Cernicharo 1993). The ratios measured in 
IRAS 18566+0408 are different from those in dark or warm 
clouds. Missing short spacing flux in the interferometer data 
can affect the observed core and envelope emissions differently, 
which in turn affects the hyperfine ratios. Therefore, we do not 
further investigate this issue quantitatively. 

3.3. Rotational Temperature 

We derive rotational temperatures of the NH3 gas. In the calcu- 
lation, we assumed LTE conditions in the gas and followed the 
procedure outlined in Ho & Townes (1983). Figure 6 presents a 
map of rotational temperature derived from the NH3 (1,1) and 

(2.2) lines. In most of the core, temperatures are around 10 K 
to 15 K. Higher rotational temperatures of about 45K are found 
toward MM-1, and along the ridge of SiO emission. In this high 
temperature region, there exists an area where the rotational tem- 
peratures cannot be derived. This is because the ratio of the (1,1) 
and (2,2) lines is sensitive to temperatures only up to 50K (Ho 
& Townes 1983). At temperatures over 50 K, the ratio of the 
two lines approaches 1 for the optically thick case, and 1.3 in 
the optically thin case for a wide range of temperatures. A small 
error in the flux measurement will result in a large uncertainty 
in rotational temperatures. Thus, the blanked area along the SiO 
outflow in Figure 6 has even higher temperatures. Assuming the 
same abundance for the ortho and para NH3 species, we use the 

(3.3) and the (1,1) lines to obtain a temperature estimate of 70 K 
for blanked area in the outflow region. 

4. Discussions 

4.1. Nature of the Continuum Peaks 

The 3 mm continuum emission reveals two peaks, MM-1 and 
MM-2, bridged by a faint extended filament. The dominant peak 
MM-1 coincides with the compact 7mm continuum source, faint 
cm continuum emission (Araya et al. 2005; Carral et al. 1999), 
and the strong peaks in NH3, SiO and HCN emission. The rota- 
tional temperature estimated from NH3 is 45 K, corresponding 
to a kinetic temperature of 80 K (Danby et al. 1988). The lumi- 
nosity of the internal source required to produce the heating can 
be estimated by the following equation (Scoville & Kwan 1976) 

T D = 65( °' lpC ) 2K4+l3 \ Lsmr ) 1 /( 1+ ' 8 >( 0,1 ) 1 /( 4+ ' 3 > J g 
r 10 5 L G / 

Here /3 is the power law index of the dust emissivity at far in- 
frared wavelengths, / = 0.08 cm 2 g _1 is the value of the dust 
emissivity at 50yum, and r the core radius. If the dust and gas 
reach a thermal equilibrium in the high density environment, 
i.e., T D = 80K at a radius r = 7000 AU, and /3 = 1, we estimate 



the luminosity of the embedded source to be 6 x 10 L©. This 
value is in rough agreement with the far IR luminosity observed 
for the region. The high luminosity, strong dust and molecular 
line emission, and high temperature in the gas all indicate em- 
bedded massive protostar(s) toward MM-1. 

The secondary mm peak MM-2 lies in a molecular ridge seen 
in NH3 and coincides with a local NH3 peak in the (1,1) and 
(2,2) transitions (see channels 85.7 and 86.4 km s -1 in Figures 4a 
and 4b). This molecular ridge appears to correspond to the dust 
emission seen at 3mm and 1.2mm. NH3 gas temperature toward 
MM-2 is about 30 K. The large amount of dense molecular gas 
and a local peak in the gas and dust emission may indicate em- 
bedded protostar(s) in the core. If the heating of the gas and dust 
is due to an internal source, we find a luminosity of 8 x 10 3 L 
for the protostar, through a similar analysis as described above 
for MM-1. On the other hand, if the heating is partially due to 
the molecular outflow in the region (see discussions below), the 
total luminosity for MM-2 would be lower. 

4.2. Massive Molecular Outflow in SiO and NH3 

The SiO emission delineates a bipolar molecular outflow in the 
region. Figure 7 shows the position-velocity plots of the SiO, 
HCN, and NH3 (2,2) and (3,3) emissions along the major axis of 
the SiO emission at a position angle of 135°. Toward the north- 
west of MM-1, the SiO emission is blue shifted with respect to 
the cloud systemic velocity of 85.2 km s _1 . The terminal veloc- 
ity of the blue-shifted SiO emission is about 15 km s _1 from the 
cloud velocity. Close to the peak MM-1, both the blue- and red- 
shifted SiO emissions are detected up to 30 km s _1 (3 cr level) 
from the cloud velocity. The blue-shifted SiO emission extends 
15" to the northwest. The red-shifted emission is far more com- 
pact spatially, with a peak detected only 2" southeast of MM-1. 

As shown in Figure 7, there exists cold and quiescent NH3 
gas along the outflow direction. The NH3 emission peaks at 
Vlsr of 85.2 km s _1 , and has a narrow FWHM of < 2 km s -1 . 
In addition to the cold gas, there exists blue-shifted high ve- 
locity emission up to V L sr of 70 km s _1 , 15 km s _1 from the 
cloud core velocity. This high velocity gas, offset to the north- 
west from the continuum peak, is part of the blue-shifted molec- 
ular outflow. The gas has an estimated temperature of 70 K. 
Although NH3 is a reliable tracer of dense gas in molecular cloud 
cores, it can be affected by molecular outflows associated with 
both low and high mass stars (LI 157: Tafalla & Bachiller 1995; 
IRAS20126+4104: Zhang et al. 1999). The high velocity NH 3 
and SiO gas has been likely accelerated and heated by shock 
processes in the outflow. 

Although the major axis of the SiO outflow agrees with that 
of the CO 2-1 outflow obtained at 11" resolution (Beuther et al. 
2002b), the CO outflow exhibits nearly symmetric bipolar mor- 
phology with the southeastern lobe much stronger than that in 
the SiO. Furthermore, the polarity of the SiO outflow appears to 
be the opposite of that of the CO: The blue-shifted SiO emis- 
sion lies to the northwest of the star, while the blue-shifted CO 
emission lies to the southeast of the star. This change of polarity 
between different tracers has been seen toward other objects (e.g. 
IRAS 20126+4104; Cesaroni et al. 1997; 1999). One plausible 
explanation is that the outflow axis lies almost in the plane of the 
sky and precesses. The low density CO gas traces the wide angle 
component in the outflow, while the high density SiO gas traces 
the well collimated jet component in the outflow. On the other 
hand, CO outflows toward massive star forming regions are of- 
ten unresolved by single dish telescopes, and break into multiple 
bipolar outflows at high angular resolution (e.g. 105358, Beuther 
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et al. 2002; AFGL 5 142, Zhang et al 2007). Furthermore, the SiO 
and CO may trace different outflows as shown in AFGL 5142 
(Zhang et al. 2007; Hunter et al. 1999) and the Orion South re- 
gion (Zapata et al. 2006). High resolution CO images of the out- 
flow will help resolving the difference. 

We compute the mass, momentum and energy in the SiO out- 
flow. Using an SiO to H2 fractional abundance of 10~ 7 (Zhang et 
al 1995), an excitation temperature of 70K derived from NH3, 
and assuming optically thin SiO emission, we obtain outflow 
mass, momentum and energy of 18 M , 200 M km s _1 , and 
4.0 x 10 46 ergs, respectively, in the blue-shifted lobe. Likewise, 
we obtain 9 M , 70 M G km s~\ and 1.0 x 10 46 ergs in the red- 
shifted lobe. Despite the uncertainty in the SiO abundance, the 
total mass, momentum and energy of 27 M Q , 270 M Q km s _1 , 
and 5.0 x 10 46 ergs are in a rough agreement (within a factor of 
2) with the estimates from the CO outflow (Beuther et al. 2002b). 

The terminal velocity of 15 km s _1 and the length of the SiO 
outflow (15") yield a dynamical time scale (Jd yn ) of 1.7 x 10 4 
years. This value is a few times smaller than that of the CO based 
on the lower angular resolution data (Beuther et al. 2002b). 
Assuming momentum conservation between the outflow and the 
underlying wind that powers the outflow: m w V w = F 'outflow, we 
can estimate the mass loss rate in the wind over the dynamical 
time scale of the outflow. The wind velocity V w can vary from 
100 km s _1 in low-mass stars to 500 km s _1 in high-mass stars 
(Zhang et al. 2005). Since the effect of inclination angle of the 
outflow is not corrected, we use a lower value of 100 km s _1 
for the wind velocity. This gives a mass loss rate (m w /Td yn ) of 
1.5 x 10~ 4 M yr _1 , and thus a lower limit to the mass accretion 
rate of 1 .5 x 10~ 4 M yr~', since some material presumably goes 
into the central protostar (Churchwell 2002). 

4.3. Heating and Weak Maser Emission 

The outflow apparently causes significant heating in the molec- 
ular gas, especially toward the blue-shifted lobe. As shown in 
Figure 6, the rotational temperature derived from the NH3 (1,1) 
and (2,2) lines is 45 K toward the position of MM-1. However, 
the temperature in the outflowing gas toward the blue-shifted 
lobe is higher, with values of 50-70 K. 

The effect of heating is further demonstrated in Figure 8. We 
compute ratios of the NH3 (3,3) and (1,1) emission. As the (3,3) 
transition has an energy level of 124 K, the higher ratios in gen- 
eral represent higher gas temperatures for thermal emission. As 
shown in Figure 8, the NH3 gas in the outflow region exhibits 
consistently higher line ratios, as expected from high tempera- 
ture regions. Toward the positions 'C and 'D', there appears to 
be heated NH3 gas in a bow shape. The NH3 spectrum toward 
'C shows red-shifted line wings, while the NH3 spectra toward 
'D' show blue-shifted line wings. Furthermore, the tips of the 
bows point away from each other, as seen in the upper-right 
panel in Figure 8. It is possible that the heating in NH3 traces 
another outflow, which is not seen in SiO. The heating likely 
arises from bow shocks as high velocity gas impinges on the 
cloud core. The potential driving source should lie in between 
'C and 'D' . However, no dust continuum emission is detected 
at a 3cr limit of 4 M . 

Figure 9 presents NH3 spectra toward positions 'A', 'B', 'C 
and 'D'. All NH3 spectra display line wing emission 15 km s _1 
blue shifted from the cloud velocity. Toward 'A', 'B' and 'C, 
the (3,3) line has a FWHM of < 1.5 km s -1 , smaller than the 
~ 3 km s" 1 FWHM in the (1,1) and (2,2) lines. The compact 
morphology, the narrow linewidth, and a lack of corresponding 



peaks in the NH3 (1,1) and (2,2) emission indicate that the com- 
pact (3,3) emission arises from population inversion, similar to 
maser emission. Maser inversion of the NH3 (3,3) has been de- 
tected toward a number of sources (e.g. W51, Zhang & Ho 1995, 
NGC6334, Kraemer & Jackson 1995, Beuther et al. 2007; DR 
21 (OH), Mangum & Wootten 1994; Mauersberger, Wilson & 
Henkel 1986; IRAS 20126+4104, Zhang et al. 2001). NH 3 (3,3) 
inversion can form through collisional excitation of NH3 by H2 
(Walmsley & Ungerechts 1983). Through collisions with H2, the 
upper level of the NH3 (3,3) (denoted as NH3 (3,3) + ) exchanges 
with its (0,0) state while the lower level of the (3,3) exchanges 
with the (1,0). Since the transition between the (3,3) + and (0,0) 
involves a change of parity and thus is more preferred, the (3,3) + 
state can be overpopulated. 

NH3 (3,3) masers are often observed in outflows. In the cases 
of IRAS 20126+4104 and NGC 6334 (Kraemer & Jackson 1995; 
Zhang et al. 2001), (3,3) masers are detected in the vicinity of 
bow shocks where outflow wind interacts with the cloud gas. 
The high velocity NH3 gas detected toward 'A', 'B' and 'C (see 
Figure 9) suggest a similar scenario. However, the spatially com- 
pact NH3 (3,3) emission appears to be resolved at an resolution 
of 1", suggesting that the emission is not strongly amplified. 

4.4. Kinematics in the NH3 Core 

In MM-1, the compact NH3 emission with a Gaussian-like pro- 
file and large linewidth distinguishes itself from the relatively 
smooth NH3 emission in the cloud core. The NH3 emission to- 
ward this position has much broader linewidth: 5.5 km s _1 at a 
spatial resolution of 3". We image the visibility data from the 
VLA C array only and obtain an angular resolution of 1". From 
this image, we obtain a fitted FWHM of 8.7 km s _1 in the (3,3) 
line. This increase indicates additional broadening in the spectral 
lines towards the inner part of the NH3 core. The NH3 compact 
structure has a size of 1.2" or 8000 AU. We estimate the mass 
in this compact structure following Ho & Townes (1983). With 
a rotational temperature of 45 K, size of 1.2" and [^-] = 10~ 7 
(Harju et al. 1993), and the assumption of LTE, we obtain a mass 
of 60 Mq. This value is consistent with the mass estimated from 
the dust emission. 

The broad NH3 linewidth toward MM-1 can arise from out- 
flow, infall/rotation or relative motion of multiple objects unre- 
solved within the synthesized beam. A collimated molecular out- 
flow is present in the SiO emission with high velocity emission 
toward the northwest and the southest of MM-1. The effect of 
the outflow also appears in the NH 3 emission, especially toward 
the northwest of MM-1 along the outflow lobe. Can the molec- 
ular outflow produce the line broadening seen in NH3 toward 
MM-1? The SiO emission is shifted from the cloud velocity. On 
the contrary, the NH3 emission peaks mostly at the cloud veloc- 
ity and appears to be Gaussian in profile. Furthermore, the mass 
and momentum in the outflow within the 5".0 x 4". 7 area (the 
synthesized beam of the SiO data) of MM-1 are 1.5 M and 10 
M km s _1 , respectively. Similarly, we compute the same quan- 
tities from the NH3 gas over a scale of 1".2, and find the mass 
and momentum of 60 M and 250 M km s _1 , respectively. The 
fractional abundances of NH3 and SiO may be uncertain and thus 
can affect the estimates provided above. Nevertheless, the com- 
parison between the masses in SiO and NH3 shows that toward 
the most central region of the core only < 3% of the material 
traced by the NH3 emission is from the molecular outflow. Since 
the outflow mass is calculated over the area (5".0 x 4" .7) 15 
times larger than that of the NH3 emission, the actual contribu- 
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tion from the outflow can be even smaller. Thus, it is unlikely 
that the molecular outflow is the main contributor to the NH3 
linewidth. 

The remaining possibilities for the large NH3 linewidth are 
motions such as infall/rotation or relative motion of multiple 
cores within the synthesized beam. Higher angular resolution 
observations in dust continuum and spectral lines will be fruit- 
ful in distinguishing these possibilities. If the NH3 linewidths 
are due to rotation and infall, similar to the signature seen in 
NH 3 toward IRAS 20126+4104 (Zhang et al. 1998), the dynam- 
ical mass, assuming gravitationally bound motion, derived us- 
ing M = is 35 Mo, for V ro , = 3 km s" 1 at R = 7000AU. 
This is compatible with the mass estimate from dust emission 
and NH3 at a similar scale. The mass infall rate, estimated using 
AnR 2 n H2 V in f a ii, is 1.5 x 10~ 3 M yr _1 for V,„/ fl « = 3 km s^ 1 and 
(is, = 10 6 cirT 3 . Thus, the mass loss rate in the wind, 1.5 x 10~ 4 
M yr~' (Section [4.2b . is 10% of the infall rate. Assuming that 
10% to 30% of the infalling mass is ejected in the outflow, the 
accretion luminosity, estimated from amounts to 3-4x 10 4 
L Q , about half of the far-IR luminosity. 

5. Conclusion 

We conducted observations of the high-mass protostellar candi- 
date IRAS 18566+0408 with the VLA and OVRO interferome- 
ters. 

(1) We resolve a collimated outflow in SiO (2-1) emission. 
The outflow is also detected in the line wings of NH3 inversion 
transitions and produces significant heating of the molecular gas 
(up to 70K). Compact features in the NH3 (3,3) line are detected 
along the outflow. The narrow linewidths of < 1 .5 km s _1 suggest 
that they are weakly amplified maser emission. 

(2) The 87 GHz emission reveals two peaks MM-1 and MM- 
2. The internal heating in the MM-1 core calls for embedded 
massive young star(s). The NH3 linewidth toward MM-1 is much 
broader than the typical linewidth of < 2 km s _1 in the extended 
core, and increases inward from 5.5 km s at the 3" scale to 
8.7 km s _1 at the 1" scale. The motion is consistent with rota- 
tion/infall, but can also arise from relative motions of unresolved 
protostellar cores. 
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Instrument 


Date of 


Line 


Bandwidth 


Spectral 


Integration 




Observations 




(MHz) 


Res. (km s~') 


Time (hr) 


VLA-CnB 


2001/07/23 


NH 3 (1,1), (2,2) 


3.12 


0.6 


0.7 


VLA-DnC 


2001/10/01 


NH 3 (3,3) 


3.12 


0.3 


1.0 


VLA-CnB 


2002/09/26 


NH 3 (1,1), (2,2) 


3.12 


0.6 


1.2 


VLA-CnB 


2002/09/26 


NH 3 (3,3) 


3.12 


0.6 


1.0 


VLA-CnB 


2002/09/26 


23 GHz 


25 




1.0 


v LA-unL, 


ZUUj/Ul/ 10 


MO /I 1 \ CO 

INrl3 (1,1 


J.VZ 




U.o 


VLA-DnC 


2003/01/24 


NH 3 (3,3) 


3.12 


0.6 


1.0 


VLA-DnC 


2003/02/04 


43 GHz 


50 




2 


OVRO-E 


2002/11/02 


HCN (1,0) 


30 


3.4 


4.0 


OVRO-E 


2002/11/02 


SiO (2-1) 


30 


1.7 


4.0 


OVRO-E 


2002/1 1/02 


HCO + (1-0) 


7.5 


0.84 


4.0 


OVRO-E 


2002/11/02 


87 GHz 


4000 




4.0 


OVRO-H 


2002/12/01 


87 GHz 


4000 




2.0 
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Fig. 1. Continuum emission at 43 GHz and 87 GHz toward IRAS 18566+0408. The contour levels are in steps of 0.25 mJy/beam 
for the 43 GHz continuum image, and 1.5 mJy/beam for the 87 GHz continuum image. The 'star' symbol and 'triangle' mark the 
continuum peaks MM-1 and MM-2, respectively. The size of the synthesized beam is marked by the shaded ellipse at the lower- left 
corner of each panel. 
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10 10 10 10 

Wavelength (/zni) 



Fig. 2. The spectral energy distribution of IRAS 18566+0408. The dashed line represents a cold dust component, the dotted line 
represents the warm dust component, the dash-dot line represents the free-free emission, and the solid line represents the sum of all 
three components. 
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Fig. 3. The integrated emission of the NH 3 (J,K)=(U), (2,2), (3,3), SiO J=(2-l), HCN J=(2-l) lines. The velocity range of the 
integration is 81 to 88 km s _1 for the NH3 (1,1) and (2,2) lines. The NH3 (3,3) and HCN lines are integrated over the entire spectral 
line. For the SiO line, the blue-shifted emission (solid contours) is integrated from 50 to 83 km s _1 and the red-shifted emission 
(dashed contours) is integrated from 88 to 1 10 km s -1 . The NH3 images are made from the VLA-D and C configuration data. The 
SiO and HCN data images are made from the OVRO E configuration data. The contour levels are in steps of 0.015 Jy km s _1 for 
the NH 3 (1,1) and (2,2) lines, 0.02 Jy km s 1 for the NH 3 (3,3), 0.3 Jy km s 1 for the SiO (2-1), and 0.6 Jy km s 1 for the HCN 
(1-0). The 'star' symbol and 'triangle' mark the continuum peaks MM-1 and MM-2, respectively. The beam sizes are plotted at the 
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Fig. 4a. The channel maps of the NH3 (J,K)=(1,1) (Fig. 4a), (2,2) (Fig. 4b) and (3,3) (Fig. 4c) lines. The images were made using 
the VLA-D configuration data only. The contour levels are in steps of 8 mJy/beam (~ 4cr) for the NH3 (1,1) and (2,2) lines and in 
steps of 9 mJy/beam (~ 3cr) for the (3,3) line. The beam sizes are plotted at the lower-left of the first panel. The velocity of the 
channel is plotted at the upper-right of each panel. The (1,1) and (2,2) images have a spectral resolution of 0.6 km s~\ while the 
(3,3)images have a spectral resolution of 0.3 km s _1 . The 'star' symbol and 'triangle' mark the continuum peaks MM-1 and MM-2, 
respectively. The crosses mark the position of the hot NH3 (3,3) features. 
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Fig. 5a. The channel maps of the SiO J=(2-l) (Fig. 5a) and HCN J=(l-0) (Fig. 5b) lines. The images were made using the OVRO E 
configuration data only. The contour levels are in steps of 0.05 Jy/beam for the SiO line, and 0.07 Jy/beam for the HCN line (~ 3<x). 
The beam sizes are plotted at the lower-left of the first panel. The velocity of the channel is plotted at the upper-right of each panel. 
The 'star' symbol and 'triangle' mark the continuum peaks MM-1 and MM-2, respectively. The crosses mark the position of the hot 
NH 3 (3,3) features. 



Zhang et al.: Jet-like Outflow toward IRAS 18566+0408 





Fig. 6. The color scale shows the rotation temperatures derived from the NH3 (1,1) and (2,2) lines. The contours represent the 
integrated emission of the SiO 2-1 line. The 'star' symbol and 'triangle' mark the continuum peaks MM-1 and MM-2, respectively. 
The crosses mark the position of the hot NH3 (3,3) features. 
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Fig. 7. The position-velocity plots of the SiO 2-1, HCN 1-0, NH 3 (J,K)=(1,1), (2,2) and (3,3) lines. The abcissa is offset from the 
position of the dust continuum with position angle (PA) of 135°. The positive offset represents the southeast of the dust continuum 
peak. The contours for the SiO and HCN are plotted in steps of +0.05 Jy/beam for the SiO line, and ±0.07 Jy/beam for the HCN 
line (~ 3o-), respectively. The NH 3 data are contoured at 3 x (1, 2, 4, 6, 8, 10, 15, 20, 25) mJy/beam. The three HCN hyperfine 
components are labeled assuming a cloud systemic velocity of 85.2 km s -1 . The NH3 (2,2) emission around 72 km s _1 is the satellite 
hyperfine component. 




Fig. 8. The ratio of the integrated NH3 emission (3,3)/(l,l) (color) overlaid on the integrated emission of the NH3 (J,K)=(1,1) line 
(contours). The velocity range of the integration in units of km s _1 are given at the upper left corner of each panel. The 'star' symbol 
and 'triangle' mark the continuum peaks MM-1 and MM-2, respectively. The crosses marks the position of hot NH3 (3,3) features. 
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Fig. 9. Spectra of the NH 3 (J,K)=(1,1), (2,2) and (3,3) transitions from the compact components A, B, C, and D. 



